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We describe the first-principles design and subsequent synthesis of a new material with the specific 
functionalities required for a solid-state-based search for the permanent electric dipole moment of 
the electron. We show computationally that perovskite-structure europium barium titanate should 
exhibit the required large and pressure-dependent ferroelectric polarization, local magnetic mo- 
ments, and absence of magnetic ordering at liquid helium temperature. Subsequent synthesis and 
characterization of Euo.sBao.sTiOa ceramics confirm the predicted desirable properties. 



The Standard Model of particle physics incorporates 
the breaking of the discrete symmetries of parity (P) 
and the combined charge conjugation and parity {CP). 
It is thought however, that the CP-violation within the 
framework of the Standard Model is insufficient to ex- 
plain the observed matter-antimatter asymmetry of the 
Universe^, therefore a so far unknown source of CP- 
violation likely exists in nature. The existence of a non- 
zero permanent electric dipole moment (EDM) of a par- 
ticle, such as an electron, neutron, or atom, would vio- 
late time reversal (T) symmetry (Fig. [TJ and therefore 
imply CP- violation through the CPT theorerrP. In the 
Standard Model these EDMs are strongly suppressed, the 
theoretical predictions lying many orders of magnitude 
below the current experimental limits. However, many 
theories beyond the Standard Model, such as supersym- 
metry, contain a number of CP- violating phases that lead 
to EDM predictions within experimental reacrP. Search- 
ing for EDMs therefore constitutes a background-free 
method of probing the CP- violating physics beyond the 




FIG. 1: Illustration that an electron with an electric dipole 
moment violates time-reversal symmetry. Both the electric 
dipole moment (+ and - symbols; orange shading) and mag- 
netic moment (blue arrow) of the electron lie along the same 
axis as the electron spin (black arrow). The operation of time 
reversal reverses the magnetic moment but does not affect the 
electric dipole moment; therefore an electron with a non-zero 
electric dipole moment violates time-reversal symmetry. 



Standard Model. 

A number of experimental EDM searches are currently 
under way or are being developed - systems studied 
in these experiments include diatomic molecules^', dia- 
magnetic atoms 6 8 , molecul ar ion ;P, cold atoms^, neu- 
trons^, liquids^, and solids^UJJ - with one of the most 
promising novel techniques being electric-field-correlated 
magnetization measurements in solidsPHUl. This tech- 
nique rests on the fact that, since spin is the only intrin- 
sic vector associated with the electron, a non-vanishing 
electron EDM is either parallel or antiparallel to its spin 
and hence its magnetic moment. As a result, when an 
electric field, which lifts the degeneracy between electrons 
with EDMs parallel and antiparallel to it, is applied to 
a sample, the associated imbalance of electron popula- 
tions generates a magnetization (Fig. [2j. The orienta- 
tion of the magnetization is reversed when the electric 
field direction is switched; in our proposed experiment 
we will monitor this change in sample magnetization us- 
ing a SQUID magnetomete r 18 ! 19 l Such magnetoelectric 
responses in materials with permanent macroscopic mag- 
netizations and polarizations are of great current interest 
in the materials science community because of their po- 
tential for enabling novel devices that tune and control 
magnetism using electric fields^. 

Since the experiment aims to detect the intrinsic mag- 
netoelectric response associated with the tiny electric 
dipole moment of the electron, the design constraints on 
the material are stringent. First, the solid must con- 
tain magnetic ions with unpaired spins, since the equal 
and opposite spins of paired electrons have corresponding 
equal and opposite EDMs and contribute no effect. Sec- 
ond, it must be engineered such that the conventional 
linear magnetoelectric tensor is zero; our approach to 
achieving this is to use a paramagnet in which the con- 
ventional effect is forbidden by time-reversal symmetry^. 
To reach the required sensitivity, a high atomic density 
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FIG. 2: Schematic of the physics underlying the experiment 
to search for the electron EDM. The energy of electrons with 
EDMs parallel to the effective electric field E* is lower than 
that for electrons with anti-parallel EDMs by an amount 
A — E* ■ d e . As a result there is a population imbalance (ex- 
aggerated for clarity in the figure), and, since the magnetic 
moments are aligned along the EDM directions, a correspond- 
ing net magnetization, M. When the electric field is reversed 
there is a magnetization reversal, AM, which can be detected 
using a sensitive magnetometer. 



of magnetic ions (n 10 22 cm -3 ) is needed, and these 
magnetic ions must reside at sites with broken inversion 
symmetry. The energy splitting A shown in Fig. [2] is 
proportional to the product of the effective electric field 
experienced by the electron, E* , and its electric dipole 
moment, d e . The effective electric field, which is equal to 
the electric field one would have to apply to a free electron 
to obtain the same energy splitting, is in turn determined 
by the displacement of the magnetic ion from the center 
of its coordination polyhedron; for a detailed derivation 
see Refp2. For example, in Euo.5Bao.5Ti0 3 ceramics (see 
below) with ~1 /iC/cm 2 remanent polarization, the mean 
displacement of the Eu 2+ ion with respect to its oxygen 
cage is 0.01 A and this results in an effective electric field 
of ~10 MV/cm, even when no external electric field is 
applied. We choose a ferroelectric so that it is possible 
to reverse the direction of the ionic displacements, and 



hence of the effective electric field, with a moderate ap- 
plied electric field. Finally, the experiment will be per- 
formed inside liquid helium, so the material properties 
described above must persist at low temperature. A de- 
tailed derivation of the dependence of the sensitivity on 
the material parameters is given in RefP^H Note that con- 
ventional impurities such as defects or domain walls are 
not detrimental to the experiment since they do not vio- 
late time-reversal symmetry. In summary, the following 
material specifications will allow a sensitive EDM search 
to be mounted: (i) The material should be ferroelectric, 
with a large electric polarization, and switchable at liquid 
He temperature, (ii) There should be a high concentra- 
tion of ions with local magnetic moments that remain 
paramagnetic at liquid He temperature; both long-range 
order and freezing into a glassy state must be avoided. 

(iii) The local environment at each magnetic ion should 
be strongly modified by the ferroelectric switching, and 

(iv) the sample should be macroscopic. With these ma- 
terial properties, and optimal SQUID noise levels, the 
projected experimental sensitivity is 10 8 e.cm after ten 
days of averaging^. 

No known materials meet all the requirements. Indeed 
the contra-indication between ferroelectricity and mag- 
netism has been studied extensively over the last decade 
in the context of multiferroics^, where the goal has been 
to achieve simultaneous ferroelectric and ferromagnetic 
ordering at high temperature. In spite of extensive ef- 
forts, a room temperature multiferroic with large and ro- 
bust ferroelectricity and magnetization at room temper- 
ature remains elusive. While the low temperature con- 
straints imposed here seem at first sight more straight- 
forward, avoiding any magnetic ordering at low temper- 
ature, while retaining a high concentration of magnetic 
ions poses a similarly demanding challenge. In addition 
the problem of ferroelectric switchability at low temper- 
ature is challenging, since coercivities tend to increase as 
temperature is lowerecPR 

We proceed by proposing a trial compound and cal- 
culating its properties using density functional theory 
to determine whether an experimental synthesis should 
be motivated. We choose an alloy of europium titanate, 
EuTi0 3 and barium titanate, BaTi03, with motivation 
as follows: To incorporate magnetism we require unfilled 
orbital manifolds of localized electrons; to avoid mag- 
netic ordering the exchange interactions should be small. 
Therefore the tightly bound 4/ electrons are likely to 
be the best choice. For conventional ferroelectricity we 
require transition metal ions with empty d orbitals to al- 
low for good hybridization with coordinating anions on 
off-centering^. (Note that while here we use a conven- 
tional ferroelectric mechanism, many alternative routes 
to ferroelectricity that are compatible with magnetism - 
and which could form a basis for future explorations - 
have been recently identified; for a review see Ref. |26| . 
Both EuTi0 3 and BaTi0 3 form in the AB0 3 perovskite 
structure, with divalent Eu 2+ or Ba 2+ on the A site, and 
formally d° Ti 4+ on the B site. BaTi0 3 is a prototypical 
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ferroelectric with a large room temperature polarization 
of 25 /iC/cm 2 .^ In the cubic paraelectric phase its lat- 
tice constant is 3.996 APh The Ba 2+ ion has an inert gas 
electron configuration and hence zero magnetic moment. 

The lattice parameter of EuTiC>3 is 3.905 A^, notably 
smaller than that of BaTiC>3. It is not ferroelectric, but 
has a large dielectric constant (e ss 400) at low tem- 
perature, indicative of proximity to a ferroelectric phase 
transition; indeed it h as re cently been reported to be 
a quantum paraelectric! 29 ! 30 !. First-principles electronic 
structure calculations have shown that ferroelectricity 
should be induced along the elongation direction by ei- 
ther compressive or tensile strairM The Eu 2+ ion has 
seven unpaired localized 4/ electrons resulting in a large 
spin magnetization of 7 fiB, and EuTiC>3 is an antifcr- 
romagnet with G-type ordering at a low Neel temper- 
ature of ~5.3Kp2Ml. (Independently of the study pre- 
sented here, EuTiC>3 is of considerable current interest 
because its dielec tric re sponse is strongly affected by the 
magnetic orderingp^Ml anc j because of its unusual third 
order magnetoelectric response^. These behaviors indi- 
cate coupling between the magnetic and dielectric orders 
caused by se nsitivi ty of the polar soft mode to the mag- 
netic orderin g 31 ! 35 !.) 

Our hypothesis is that by alloying Ba on the A-site 
of EuTi03, the magnetic ordering temperature will be 
suppressed through dilution, and the tendency to fer- 
roelectricity will be increased through the expansion of 
the lattice constant. Our hope is to identify an alloying 
range in which the magnetic ordering temperature is suf- 
ficiently low while the ferroelectric polarization and the 
concentration of magnetic ions remain sufficiently large. 
In addition, we expect that the polarization will be sen- 
sitive to the lattice constant, allowing its magnitude and 
consequently the coercivity, to be reduced with pressure. 



I. FIRST-PRINCIPLES CALCULATIONS 

Taking the 50/50 (Eu,Ba)Ti03 ordered alloy as our 
starting point (Fig. [3] inset), we next calculate its proper- 
ties using first-principles. For details of the computations 
see the Methods section. 

We began by calculating the phonon dispersion for the 
high symmetry, cubic perovskite reference structure at a 
lattice constant of 3.95 A (chosen, somewhat arbitrarily, 
for this first step because it is the average of the ex- 
perimental BaTiOs and EuTi03 lattice constants), with 
the magnetic spins aligned ferromagnetically; our results 
are shown in Fig. |3j plotted along the high symmetry 
lines of the Brillouin zone. Importantly we find a po- 
lar T-point instability with an imaginary frequency of 
103i cm -1 which is dominated by relative oxygen - Ti/Eu 
displacements (the eigenmode displacements for Eu, Ba, 
Ti, 0|| and ± are 0.234, -0.059, 0.394, -0.360 and -0.303 
respectively); such polar instabilities are indicative of a 
tendency to ferroelectricity. The zone boundary rota- 
tional instabilities that often occur in perovskite oxides 
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FIG. 3: Calculated phonon dispersion of ferromagnetic 
Euo.sBao.sTiOs in its high symmetry reference structure with 
pseudo-cubic lattice constant ao = 3.95 A. The imaginary fre- 
quency polar phonon at T indicates a structural instability to 
a ferroelectric phase. The inset shows the supercell of the 
ferromagnetic Euo.sBao.sTiOa ordered alloy used in our cal- 
culations. The Ti and O ions are omitted for clarity; arrows 
represent the Eu magnetic moments. 



and lead to non-polar, antiferrodistortive ground states 
are notably absent (in fact the flat bands at ^60 cm -1 
are stable rotational vibrations). Interestingly we find 
that the Eu ions have a significant amplitude in the soft- 
mode eigenvector, in contrast to the Ba ions both here 
and in the parent BaTiOs. 

Next we performed a structural optimization of both 
the unit cell shape and the ionic positions of our 
Euo.sBao.sTiOs alloy with the total volume constrained 
to that of the ideal cubic structure studied above (3.95 3 
A 3 per formula unit). Our main finding is that the 
Euo.sBao.sTiOs alloy is polar with large relative displace- 
ments of oxygen and both Ti and Eu relative to the high 
symmetry reference structure. Using the Berry phase 
method we obtain a ferroelectric polarization value of 
P = 23 /iC/cm 2 . Our calculated ground state is or- 
thorhombic with the polarization oriented along a [011] 
direction and lattice parameters a — 3.94 A, b = 5.60 
A and c = 5.59 A. As expected from our analysis of the 
soft mode, the calculated ground state is characterized by 
large oxygen - Ti/Eu displacements, and the absence of 
rotations or tilts of the oxygen octahedra. Importantly, 
the large Eu amplitude in the soft mode manifests as a 
large off-centering of the Eu from the center of its oxy- 
gen coordination polyhedron in the ground state struc- 
ture. The origin of the large Eu displacement lies in its 
small ionic radius compared with that of divalent Ba 2+ : 
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The large coordination cage around the Eu ion which is 
imposed by the large lattice constant of the alloy results 
in under-bonding of the Eu that can be relieved by off- 
centering. Indeed, we find that in calculations for fully 
relaxed single phase EuTiC>3, the oxygen octahedra tilt 
to reduce the volume of the A site in a similar manner to 
those known to occur in SrTi03, in which the A cation 
size is almost identical. This Eu off-centering is desirable 
for the EDM experiment because the change in local en- 
vironment at the magnetic ions on ferroelectric switching 
determines the sensitivity of the EDM measurement. 



volume (A 3 ) 


P (AtC/cm 2 ) 


61.63 (constrained) 


23 


62.30 (experimental) 


28 


64.63 (relaxed) 


44 



TABLE I: Calculated ferroelectric polarizations, P, of 
Euo.sBao.sTiOs at three different volumes. 

We note that the magnitude of the polarization is 
strongly dependent on the volume used in the calcula- 
tion (Table [I]). At the experimental volume (reported in 
the next section), which is only slightly larger than our 
constrained volume of 3.95 3 A 3 , we obtain a polariza- 
tion of 28 /LtC/cm 2 . At full relaxation, where we hnd a 
larger volume close to that of BaTiC>3, we obtain a po- 
larization of 44 /iC/cm 2 , almost certainly a substantial 
over-estimate. This volume dependence suggests that the 
use of pressure to reduce the lattice parameters and sup- 
press the ferroelectric polarization could be a viable tool 
for reducing the coercivity at low temperatures. Indeed 
our computations show that, at a pressure corresponding 
to 2.8 GPa applied to the experimental volume the theo- 
retical structure is cubic, with both the polarization and 
coercive field reduced to zero. 

Finally, to investigate the likelihood of magnetic order- 
ing, we calculated the relative energies of the ferromag- 
netic state discussed above and of two antifcrromagnetic 
arrangements: planes of ferromagnetically ordered spins 
coupled antiferromagnetically along either the pseudo- 
cubic z axis or the x or y axes. (Note that these are de- 
generate in the high-symmetry cubic structure). For each 
magnetic arrangement we re-relaxed the lattice parame- 
ters and atomic positions. As expected for the highly lo- 
calized Eu 4/ electrons on their diluted sublattice, the en- 
ergy differences between the different configurations are 
small - around 1 meV per 40 atom supercell - suggesting 
an absence of magnetic ordering down to low tempera- 
tures. While our calculations find the ferromagnetic state 
to be the lowest energy, this is likely a consequence of our 
A-site ordering and should not lead us to anticipate fer- 
romagnetism at low temperature (Note that, after com- 
pleting our study, we found a report of an early effort 
to synthesize (Eu,Ba)TiOgp^in which a large magnetiza- 
tion, attributed to A-site ordering and ferromagnetism, 
was reported. A-site ordering is now known to be dif- 
ficult to achieve in perovskite-structure oxides, however, 



and we find no evidence of it in our samples. Moreover 
the earlier work determined a tetragonal crystal structure 
in contrast to our refined orthorhombic structure.) 

In summary, our predicted properties of the 
(Eu,Ba)Ti03 alloy - large ferroelectric polarization, re- 
ducible with pressure, with large Eu displacements, and 
strongly suppressed magnetic ordering - meet the criteria 
for the electron electric dipole moment search and moti- 
vate the synthesis and characterization of the compound, 
described next. 



II. SYNTHESIS 

Euo^Bao^TiOs was synthesized by solid-state reaction 
using mechanochemical activation before calcination. For 
details see the Methods section. The density of the sin- 
tered pellets was 86-88% of the theoretical density. X-ray 
diffraction at room temperature revealed the cubic per- 
ovskite Pm3m structure with a=3. 9642(1) A. At 100 K 
we obtain an orthorhombic ground state with space group 
Amm2, in agreement with the GGA+J7 prediction, and 
lattice parameters 3.9563(1), 5.6069(2) and 5.5998(2) A. 

III. CHARACTERIZATION 

The final step in our study is the characterization of 
the samples, to check that the measured properties are 
indeed the same as those that we predicted and desired. 
Figure [4] shows the temperature dependence of the com- 
plex permittivity between 1 Hz and 1 MHz, measured us- 
ing an impedance analyzer ALPHA- AN (Novocontrol) . 
The low-frequency data below 100 kHz are affected above 
150 K by a small defect-induced conductivity and re- 
lated Maxwell- Wagner polarization; the high-frequency 
data clearly show a maximum in the permittivity near 
T C =213K indicating the ferroelectric phase transition. 
Two regions of dielectric dispersion - near 100 K and be- 
low 75 K - are seen in tan^(T); these could originate from 
oxygen defects or from ferroelectric domain wall motion. 

Measurement of the polarization was adversely affected 
by the sample conductivity above 150 K, but at lower 
temperatures good quality ferroelectric hysteresis loops 
were obtained (Fig.|4j inset). At 135 K we obtain a sat- 
uration polarization of ~8 /LtC/cm 2 . The deviation from 
the predicted value could be the result of incomplete sat- 
uration as well as the strong volume dependence of the 
polarization combined with the well-known inaccuracies 
in GGA+C7 volumes. As expected, at lower temperatures 
the coercive field strongly increases, and only partial po- 
larization switching was possible even with an applied 
electric field of 18kV/cm (at higher electric field dielec- 
tric breakdown was imminent). The partial switching is 
responsible for the apparent decrease in saturation po- 
larization below 40 K. 

Time-domain THz transmission and infrared reflectiv- 
ity spectra (not shown here) reveal a softening of the 
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FIG. 4: Temperature dependence of permittivity and dielec- 
tric loss in Euo.sBao.sTiOs ceramics. The arrows indicate the 
direction of increasing frequency and the colors are for clarity 
to assist the eye in distinguishing the lines. The inset shows 
ferroelectric hysteresis loops measured at three temperatures 
and 50 Hz. 



polar phonon from ^40 cm" 1 at 300 K to ^15 cm -1 at 
T Cl and then its splitting into two components in the 
ferroelectric phase. Both components harden on cooling 
below T c , with the lower frequency component remaining 
below 20 cm" 1 down to 10 K, and the higher-frequency 
branch saturating near 70 cm -1 at 10 K. This behavior 
is reminiscent of the soft-mode behavior in BaTiOj^l 
However, when we extract the contribution to the static 
permittivity that comes from the polar phonon, we find 
that it is considerably smaller than our measured value 
(Fig. [4j indicating an additional contribution to the di- 
electric relaxation. Our observations suggest that the 
phase transition is primarily soft-mode driven, but also 
exhibits some order-disorder character. 

Finally, we measured the magnetic susceptibility \ at 
various static magnetic fields as a function of temper- 
ature down to 0.4 K. (For details see the Methods sec- 
tion.) Our results are shown in Fig. [5] x(T) peaks 
at T ~1.9K indicating an absence of magnetic order- 
ing above this temperature. The x(^) data up to 300 K 
show Curie- Weiss behavior x(T) = jf^g with 0=-1.63 K 
and C — 0.017 emuK/(gOe). The peak in susceptibil- 
ity at 1.9K is frequency independent and not influenced 
by zero field heating measurements after field cooling, 
confirming antiferromagnetic order below T/v = 1.9 K. 
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FIG. 5: Temperature dependence of ac magnetic suscepti- 
bility, x> a t various static magnetic fields and frequency of 
214 Hz. Inset shows magnetization curves at various tem- 
peratures. We note that no hysteresis in magnetization was 
observed. 



As in pure EuTi0 3 , the x(T) peak is suppressed by a 
static external magnetic field, indicating stabilization of 
the paramagnetic phased. Magnetization curves (Fig. [5] 
inset) show saturation above 2 x 10 4 Oe at temperatures 
below Tjq and slower saturation at 5 K. No open magnetic 
hysteresis loops were observed. 

In summary, we have designed a new material - 
Euo.sBao.sTiOs - with the properties required to en- 
able a measurement of the EDM to a higher accuracy 
than can currently be realized. Subsequent synthesis of 
Euo.sBao.sTiOs ceramics confirmed their desirable fer- 
roelectric polarization and absence of magnetic ordering 
above 1.9 K. The search for the permanent dipole mo- 
ment of the electron using Euo.5Bao.5Ti03 is now un- 
derway. Initial measurements have already achieved an 
EDM upper limit of 5 xlO -23 e.cm, which is within a 
factor of 10 of the current record with a solid-state-based 
EDM searctM We are currently studying a number of 
systematic effects that may mask the EDM signal. The 
primary error originates from ferroelectric hysteresis- 
induced heating of the samples during polarization rever- 
sal. This heating gives rise to a change in magnetic sus- 
ceptibility, which, in a non-zero external magnetic field, 
leads to an undesirable sample magnetization response. 
We are working to control the absolute magnetic field 
at the location of the samples to the 0.1 /LtG level. Our 
projected sensitivity of 10~ 28 e.cm should then be achiev- 
able. 
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VI. METHODS 

A. Computational details 

We performed first-principles density-functional calcu- 
lations within the spin-polarized generalized gradient ap- 
proximation 

(GGAJpH. The strong on-site correlations 
of the Eu 4/ electrons were treated using the GGA+/7 
methocP^ with the double counting treated within the 
Dudarev approactP^ and parameters U = 5.7 eV and 
J = 1.0 eV. For structural relaxation and lattice dy- 
namics we used the Vienna Ab Initio Simulation Pack- 
age (VASPj^with the default projector augmented-wave 
(PAW) potentials^ (valence-electron configurations Eu: 



5s 2 5p 6 4/ 7 6s 2 , Ba: 5s 2 5p 6 6s 2 , Ti: 3s 2 3p 6 3d 2 4s 2 and O: 
2s 2 2p 4 .) Spin-orbit interaction was not included. 

The 50/50 (Eu,Ba)Ti0 3 alloy was represented by an 
ordered A-site structure with the Eu and Ba ions alter- 
nating in a checkerboard pattern (Fig. |3j inset). Struc- 
tural relaxations and total energy calculations were per- 
formed for a 40-atom supercell (consisting of two 5-atom 
perovskite unit cells in each cartesian direction) using 
a 4 x 4 x 4 T-centered fc-point mesh and a plane-wave 
cutoff of 500 eV. Ferroelectric polarizations and Born 
effective charges were calculated using the Berry phase 
method^. Lattice instabilities were investigated in the 
frozen-phonon scheme^ 4 ^ for an 80 atom supercell using 
a T-centered 2x2x2 /c-point mesh and 0.0056 A atomic 
displacements to extract the Hellman-Feynman forces. 
B. Synthesis 

EU2O3, Ti02 (anatase) and BaTiOs powders (all from 
Sigma-Aldrich) were mixed in stoichiometric ratio then 
milled intensively in a planetary ball micro mill Fritsch 
Pulverisette 7 for 120 min. in a dry environment followed 
by 20 min. in suspension with n-heptane. Zr02 grinding 
bowls (25 ml) and balls (12 mm diameter, acceleration 
14 g) were used. The suspension was dried under an IR 
lamp and the dried powder was pressed in a uniaxial press 
(330 MPa, 3 min.) into 13 mm diameter pellets. The pel- 
lets were calcined in pure H 2 atmosphere at 1200 °C for 
24 hr (to reduce Eu 3+ to Eu 2+ ), then milled and pressed 
by the same procedure as above and sintered at 1300 °C 
for 24 hr in Ar + 10% H2 atmosphere. Note that pure H2 
can not be used for sintering without adversely increasing 
the conductivity of the sample. 



C. Characterization 

Magnetic susceptibility was measured using a Quan- 
tum Design PPMS9 and a He 3 insert equipped with a 
home-made induction coil that allows measurement of ac 
magnetic susceptibility, \ from 0.1 to 214 Hz. 
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